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Abstract 



We present the first numerical simulation of actin-driven propulsion by 
elastic filaments. Specifically, we use a Brownian dynamics formulation of 
the dendritic nucleation model of actin-driven propulsion. We show that 
the model leads to a self- assembled network that exerts forces on a disk 
and pushes it with an average speed. This simulation approach is the first 
to observe a speed that varies non-monotonically with the concentration 
of branching proteins (Arp2/3), capping protein and depolymerization rate 
(ADF), in accord with experimental observations. Our results suggest a 
new interpretation of the origin of motility that can be tested readily by 
experiment. 
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Introduction 

There is a type of biological motility, used in a form of cell crawling (1) and 
by intracellular pathogens such as Listeria monocytogenes(2) , that is driven 
not by motor proteins but by biological self-assembly of the protein actin. 
During this process, ATP hydrolysis and activation of the protein complex 
Arp2/3 drive actin self-assembly from monomers (G-actin) to branched net- 
works of filaments (F-actin) (3), thus providing the necessary thermody- 
namic free energy to push a bacterium or a cell forward (4,5). This driven, 
non-equilibrium self-assembly process is regulated by a cadre of proteins. 
It is now possible to drive a latex bead through a buffer solution contain- 
ing only these proteins (6-9). Such beads travel through solution propelled 
by a dense branched actin network at their rear, demonstrating that non- 
equilibrium self-assembly of F-actin is sufficient to drive motility. 

The standard biochemical model for the regulation of actin-self-assembly- 
driven motility is the dendritic nucleation model (3,10,11). In this model, 
actin self-assembles (or "polymerizes") into filaments preferentially at one 
end (the barbed end) and " depolymerizes" preferentially at the other end 
(the pointed end). Proteins such as WASP at the moving surface (the rear 
end of the Listeria bacterium or moving latex bead, or the leading edge of 
the membrane of a crawling cell) recruit and activate the Arp2/3 protein 
complex. The activated Arp2/3 catalyzes the nucleation of new branches 
from pre-existing actin filaments, thus creating new growing barbed ends 
near the moving surface. To sustain motion, two other essential proteins 
regulate the turnover of actin monomers: severing protein (ADF), which 
raises the depolymerization rate by severing filaments in two, and capping 
protein (Cap), which covers barbed ends and prevents further growth. Thus, 
filaments just behind the moving surface at the front of the branched net- 
work tend to grow due to Arp2/3 and WASP, while filaments at the far end 
of the branched network tend to depolymerize away due to ADF and Cap. 

A key physical question that arises is: how does the self-assembly of a 
branched network generate forces and produce motion? Many models have 
been developed to show how the polymerization of a single actin filament can 
produce a force (12-22). Other models show how the dendritic nucleation 
model creates a branched network morphology (23-28). Relatively few mod- 
els have considered how polymerization of a branched network might lead to 
force generation; of these, some treat the network as an elastic continuum 
(29-31). Only three approaches explicitly incorporate the morphology of the 
dendritic nucleation model to produce force and motion (28,32,33). In all 
three of these simulation models, mass is not conserved; monomers spring 
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into existence and become capable of exerting forces only when they join 
filaments, and vanish when they fall off. As a result, matter is created just 
behind the moving surface, leading to motility as an artifact. 

In this paper, we use Brownian dynamics to demonstrate that force and 
motion can indeed emerge from the growth of a branched network in a 
physically-consistent model. We demonstrate that our model is the first 
to capture key properties of the dendritic nucleation model by reproducing 
the characteristic dependence of speed on the concentrations of Arp2/3, 
capping protein, ADF and actin (8,34,35). Our simulation suggests a new 
understanding of the mechanism driving motility: the disk emits activated 
Arp2/3 complex, which gives rise to a buildup of F-actin just behind the 
disk. If there is a repulsive interaction between the disk and the actin, the 
disk will move forward to avoid the actin recruited by Arp2/3. We propose 
explicit experiments to test this new picture. 

Simulation Model 

Relevant time scales within the dendritic nucleation model span six orders 
of magnitude. The longest time scale is set by kinetic events such as the 
depolymerization rate (~ Is) (3), while the shortest time scale is deter- 
mined by diffusion and collision of G-actin monomers (~ 1/xs) and the high 
frequency dynamics of filaments. The wide range of important time scales 
poses a challenge to computer simulation. Previous approaches avoid this 
problem by treating free monomers and those in filaments very differently, 
leading to potential artifacts (36). If one insists on treating free monomers 
and monomers within filaments consistently, one must use a time step that 
is small enough to capture their short-time-scale dynamics when integrat- 
ing the equations of motion. On the other hand, in order to study the 
steady-state, one must be able to reach time scales that are long compared 
to the slowest reaction rate involved. The compromise that we have chosen 
is to narrow the range of time scales by increasing the slowest rates, such 
as the depolymerization rate, by 5 orders of magnitude and and decreasing 
the filament stiffness by one order of magnitude (see Table 1). We also use 
the enhanced depolymerization rate to mimic the action of severing protein 
(ADF). The details of our model are presented in the Methods section. 

To offset some of these changes, we adjust other variables so that the 
steady state fluxes are comparable to those observed experimentally. For 
example, to offset the effect of our artificially-high depolymerization rate, 
we increase the typical concentration of the G-actin monomers such that 
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the ratio of the effective polymerization rate and the depolymerization rate, 
.K+[G-actin]/ K_ , is close to the typical experimental value. 

In testing our model, our aim is not to reproduce numerically accurate 
results but to capture experimentally-observed trends and understand what 
factors control them. In particular, our goal is to gain insight into the 
mechanism that leads to motility. We will show that the origin of motility 
in our system (see Results) suggests a possible mechanism for the real system 
that yields a reasonable speed within a simple order-of-magnitude estimate 
(see Discussion). 

In order to check that our conclusions do not result from the unphysical 
parameters we have chosen, we have varied the parameters over a range. 
For example, most of our runs were carried out for a bending stiffness of 
Kb = WOk^T (see Eq. [4]), corresponding to a persistence length of 0.5 fim. 
However, we have also shown that when all other parameters are held fixed, 
we obtain the same speed for Kb = 1000A;bT, corresponding to a more 
realistic persistence length of 5/im. 

We have also checked the dependence of our results on K_ and other 
slow rates by decreasing them and showing that the trends remain the same. 

Bulk system 

Our simulation model is described in the Methods section. For systems 
that are spatially isotropic on average, we have shown that the Brownian 
dynamics results for morphology are in quantitative agreement with a mean- 
field formulation of the dendritic nucleation model (27). This mean- field 
formulation was in turn shown to be in quantitative agreement with in vitro 
experiments (37). Thus, our model yields reasonable results for the steady- 
state bulk system. 

Motility 

We now break symmetry by introducing a moving surface in the form of a 
disk, whose back surface (facing the — z direction) emits Arp2/3. This drives 
self-assembly of a branched network behind the disk, which pushes the disk 
in the +z direction. We typically begin each run with 5-10% of the actin 
monomers in dimer form and the rest as free monomers. We begin with some 
dimers as protofilaments because spontaneous nucleation of filaments, which 
occurs at a very low rate experimentally (3), is not allowed in our model. We 
find that the results are not sensitive to the fraction of initial dimers. The 
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dimers and free monomers are initially distributed randomly in the system. 
Fig. Qla) shows the displacement, z, of the disk as a function of time for a 
typical simulation run for a filament stiffness Kf, = 100 (solid curve). The 
dashed vertical lines mark the times corresponding to the snapshots (I-III). 
In the snapshots, free monomers are not shown. The black box corresponds 
to the simulation box; we have shown part of the periodic images to the 
right and left. Snapshot I displays the system at t = 70/is. At this time, 
the disk is still very close to its starting position. The dimers have grown 
into short filaments and are dispersed throughout the box. By t = 700/zms 
(snapshot II), a branched F-actin network has formed behind the disk and 
the disk has moved slightly. By t = 2100ms (snapshot III), the disk has 
moved to the right by nearly a third of the simulation box. 

Fig.[T]Ja) shows that once the disk starts moving, the trajectory is linear. 
Because there are significant fluctuations in the displacement (9,38), we 
extract speeds from trajectories that are at least 7000//S long (several times 
longer than that shown in Fig. [1]), and average over the final 3500-4200//S of 
the trajectory (it takes roughly 1000/US to reach steady state). The error bars 
for the speed in all of our figures were obtained from the standard deviation 
calculated over 5 separate simulations run under standard conditions (see 
Table 1). 

The typical speed for our simulated systems is 60/xm/s. Our speed is 
simply determined by the polymerization rate. We use K + = 63/iM _1 s~ 1 
(see Table 1) . To convert this to a net polymerization speed v p , we must 
multiply by the monomer size, a = 5nm, by the free monomer concentration 
just behind the disk, and a factor that characterizes the structure of the 
network. A reasonable approximation to this factor is cos 6, where 6 is the 
angle between the average tangent vector of filaments just behind the disk 
and the normal to the disk (13). For the conditions corresponding to Fig. [IJ 
[G-actin] ~1 mM and cos((9) = 0.1 — 0.2. This yields an estimate of the 
polymerization speed v p ~ 31 — 63/im/s, in good agreement with our result. 

The speed found experimentally is significantly slower, with a typical 
value of a fraction of a micron per minute (7-9). We find that when we 
decrease the depolymerization rate and G-actin concentration by a factor 
of 10, leaving the ratio if+[G — actin]/iT_ fixed, the speed decreases by a 
factor of ~ 10. As Table I shows, the value of K + [G — &ctm]/K^ that we 
use is close to the experimental value, but K_ and [G-actin] are much higher 
in our simulation. We would therefore expect our speed to be too high. 

It is also possible that part of the difference between our simulated speed 
and the experimental speed may be due to our neglect of filament bind- 
ing to the moving surface. Experimentally, it is known that filaments in 



the branched network bind to the proteins on the disk that active Arp2/3 
complex (14, 29, 38-41). Finite element simulations (42) suggest that the 
inclusion of a binding energy between filaments and the disk slows down the 
speed significantly and enhances fluctuations around the average speed. 

An important and surprising result of our calculation is that the speed 
is independent of bending stiffness. This is shown in Fig. [2(a), where the 
speed is the same for systems with bending stiffnesses of Kb = 1000/cb7\ 
Kb = IOO/^bT and Kb = Ok^T. Note that in the flexible case, we have not 
shown the initial start-up of the disk, which is substantially longer than for 
stiffer filaments. Our results show that flexible and stiff filament networks 
exert comparable forces as the filaments polymerize. In all cases, the speed 
is simply the polymerization speed v p . Thus, the physical origin of motility 
does not depend sensitively on the bending stiffness of the filaments, as is 
commonly believed (see Discussion). 

Origin of Steady State Force 

When the system is in steady state and there is a net force on the disk moving 
it forwards, there must be an equal and opposite net force on the actin. We 
have calculated the average force density F z (z) in the z-direction at different 
distances z from the disk (recall that the disk is constrained to move only in 
the z-direction). Fig. [2] shows F z (z) in the frame of the moving disk, where 
z = (marked by a vertical dashed line) always marks the position of the 
disk. Just behind the disk at z < 0, the force on monomers (free or bound 
in filaments) is large and negative, as expected because the force exerted 
by these monomers on the disk is positive. Note that this negative force 
persists out to about 50nm behind the disk before it drops nearly to zero. 
For z > 0, the force near the disk is positive because monomers immediately 
in front are pushed along by the disk. We have verified that the total average 
force, J F z (z)dV, exerted on the actin is equal and opposite to the force on 
the disk, as it must be. The force on the bead is on the order of 0.1 pN. 
Although this force seems small, we note that it is the magnitude of the 
force required to push a lfim bead at the experimentally-observed speed, 
and is also, by construction, the force needed to push the disk at the speed 
that we observe for the viscosity chosen; we have confirmed that the average 
force on the disk is related to its average speed by the drag on the disk, £d> 
as expected. 

Note that while the negative force extends to 50nm behind the disk, 
the total length of the actin comet tail in our simulations is about 150nm 
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(Fig. [TJIII). Thus, only a relatively small fraction of the network directly 
behind the disk is subjected to a significant backwards force. This result is 
consistent with the experimental finding that the actin network in the tail 
is stationary (7, 43). 

We remark that the force profile shown in Fig. [2](a) does not contradict 
the experimental observation that the shape of the tail can be deformed at 
distances far greater than 50nm from the surface (31), because the moving 
surface was curved in the experiment and the tail expanded as it moved 
backwards away from the surface, due to entropy or elastic stresses. 

The force on the disk can be viewed as the Newton's third law reaction 
force to the force in Fig. [2(a) on the actin, evaluated at the surface of 
the disk. Thus, uncovering the origin of the force profile behind the disk 
should help us to understand motility. The solid curve in Fig. [2](b) shows 
the density profile p(z) of actin (note that the free monomer density is 
nearly constant, with a small dip just behind the disk, so that most of the 
variation is due to monomers in filament form). In equilibrium, similar 
density profiles can arise from attraction to the surface. In that case, the 
chemical potential must be the same everywhere. However, in this steady- 
state driven system, the density profile does not arise from attractions-the 
interaction of actin with the disk is purely repulsive. Rather, the density 
profile is a non- equilibrium effect, arising from the action of Arp2/3, which 
is emitted from the disk. (In the real system, Arp2/3 is activated at the 
surface of the disk, so the disk serves as a source of activated Arp2/3.) 
The non-equilibrium density profile leads to a pressure gradient, dp/dz = 
(dp /dp) (dp /dz) = — (1/ np)dp/dz, where k is the local compressibility of the 
branched network. The importance of the compression modulus has been 
emphasized in previous models (29-31). In our case, the force generated 
depends not only on k but on the concentration gradient, dp/dz. Note that 
the pressure gradient is equal and opposite to the force per unit volume on 
the actin, shown in Fig. |2)(a). The vanishing of the force near z = — 30nm 
therefore corresponds to the maximum in the concentration there, where 

dp/dz = (Fig. ED (b)). 

We emphasize that this is not a simple osmotic pressure effect due to 
free monomers. The density of free monomers (dotted curve in Fig. [2](b)) 
is nearly constant, so that the density gradient arises from F-actin, not 
G-actin. 

The fact that the speed corresponds to the polymerization speed for 
different filament stiffness suggests that the system adjusts the force exerted 
on the disk to maintain the speed at the polymerization speed, at least at 
the small loads studied here. Therefore, it should be possible to understand 
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the underlying mechanism for motility without involving the force. The 
following interpretation does not invoke forces explicitly, but is equivalent 
to the above arguments and much simpler. In this picture, Arp2/3 complex 
recruits F-actin to the vicinity of the disk. The interaction between actin and 
the disk is repulsive, so the disk moves forward to lower the concentration 
of actin near the surface. This leads to the steady-state density profile of 
Fig. EJb) as well as steady-state motion of the disk at the net polymerization 
speed. 

Dependence on protein concentrations 

One key observation of experiments is that the speed is a non-monotonic 
function of the concentrations of the regulatory proteins involved in the 
dendritic nucleation model, namely Arp2/3, capping protein and severing 
protein. Our simulation model is the first to capture this behavior and to 
explain the physical origin of the non-monotonicity. 

Fig. EJa) shows that the speed is a non-monotonic function of Arp2/3 
concentration. Similar non-monotonic behavior has been found experimen- 
tally (8). The behavior can be understood as follows. At high Arp2/3 
concentrations, most of the excess Arp2/3 is trapped in filaments in the 
network, forming stubby branches. These short, stubby branches do little 
to increase the actin concentration behind the disk. However, they do repel 
actin monomers, lowering the concentration of free monomers at the surface 
so that fewer of them are available for polymerization. This crowding effect 
is captured for the first time in our simulation because we treat monomers 
explicitly. At high Arp2/3 concentration, there appear to be two effects 
that reduce the speed: first, the maximum in the density profile in Fig. EJ^b) 
broadens as stubby branches proliferate. Second, the concentration of G- 
actin at the surface decreases. With increasing Arp2/3 concentration, the 
G-actin concentration at the surface drops below its critical value for poly- 
merization and/or the density gradient in F-actin vanishes; at this point, 
the speed drops to zero. 

The open symbols in Fig. [3][a) show the speed as a function of [Arp2/3] 
in the case where the depolymerization rate, debranching rate and G-actin 
concentration have all been decreased by a factor of 10 and the capping 
rate has been decreased by a factor of 5, relative to the values in Table 
1. While still high, the difference between the closed and open symbols 
shows the trend to be expected if we could reduce the parameters to their 
experimental values. The overall trends are the same in both cases, but 
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the speed is slower, as discussed earlier, and the maximum speed is at lower 
Arp2/3 concentration, as one might expect. The maximum is much narrower 
as a function of [Arp2/3], which is more consistent with experimental results 
(8). 

We note that for the typical parameters listed in Table 1 as well as the 
reduced parameters, the net polymerization rate is comparable to that in the 
real system. As a result, the filament density in the comet tail is comparable 
to that observed experimentally. We find a filament density of order a few 
mM in the comet tail for our standard runs, and of order 0.35 mM for the 
runs with the reduced reaction rates. This shows that there is some decrease 
in the amount of F-actin in the comet tail with decreasing depolymerization 
rate, but the values we find compare reasonably well with previous results of 
Carlsson (28). One can also estimate the filament density from the Young's 
modulus, measured to be Y = 10 3 Pa (29). The Young's modulus for a 
network of semiflexible polymers with persistence length £ p and mesh size 
U is (44, 45) 

y = k B Tiye m (i) 

where £ m = l/\fac, where a is the filament diameter and c is the monomer 
concentration. This yields c ~ ImM, as well. 

Since the filament density in our simulation is approximately the same 
as that in experiments, it is reasonable that the monomer concentration 
should be reduced near the surface relative to its value in the bulk in the real 
system. This reduction inevitably leads to a reduction of the polymerization 
rate with increasing Arp2/3 concentration. 

Fig. [3](b) shows that the speed is also non-monotonic with capping pro- 
tein, in agreement with experiment (8, 34). In this case, it is obvious that 
too much capping will lead to a vanishing speed. If the capping rate is too 
low, however, the speed also vanishes. This is because capping and branch- 
ing act synergistically. Capping stops free barbed ends from growing, thus 
forcing the system to favor branching to generate new growing ends instead 
of merely lengthening existing filaments (46, 34, 27). The capping rate at 
the maximum of the curve in Fig. E^b) is comparable to the debranching 
rate. 

Fig. 0(c) shows the dependence of speed on the depolymerization rate. 
The closed circles correspond to the case in which the Arp2/3 protects the 
pointed end from depolymerization once it reaches a branch point, and 
prevents the branch from falling off. There is experimental evidence that 
Arp2/3 protects the pointed end from depolymerization (10, 47). In this 
case, Fig. [3jc) shows that the speed saturates with increasing depolymeriza- 
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tion rate. The open circles correspond to the case in which depolymerization 
proceeds through the branch point, and the branch falls off. This is consis- 
tent with experiments that show that when ADF cofilin is present, Arp2/3 
no longer protects the pointed end from depolymerization (47). Fig. [3|c) 
shows that in this case, the speed is non-monotonic and decreases with suf- 
ficiently high depolymerization rate. The experiments of Loisel et al. (8) 
exhibit non-monotonic dependence, similar to the open circles in Fig. [3|c). 
This suggests that ADF does indeed prevent Arp2/3 from protecting the 
pointed end from depolymerization. Note that the two curves are the same 
at low K_, and begin to deviate from each other near the maximum. This 
corresponds to where K~ is comparable to the debranching rate. 

Finally, Fig. [3](d) shows the dependence on the overall actin concentra- 
tion. Again, our results are qualitatively consistent with those of experiment 
(35). The speed increases with [Actin], because the polymerization rate in- 
creases, and saturates at high [Actin] at a maximum polymerization speed, 
v p . At high [Actin], the concentration of free actin monomers at the sur- 
face, needed for polymerization, saturates at roughly 1 mM. This saturation 
apparently occurs because the branched network becomes denser and more 
difficult for the free monomers to penetrate in order to reach the disk (24,48, 
49). 

Discussion 

Our simulations show that a physically reasonable formulation of the den- 
dritic nucleation model can lead to motility. We have taken great care to 
avoid possible artifacts. For example, we treat free monomers at the same 
level as monomers in filaments so that there is no artificial mass transfer or 
dynamical discontinuity when monomers join or leave filaments. 

We have shown that the speed does not depend on the bending stiffness 
of the filaments (33, 50). This surprising result appears to be consistent 
with the observation that amoeboid sperm of nematodes (50, 51) moves 
using a structurally different filament composed of major sperm protein 
(MSP) instead of actin. These MSP filaments assemble into thick bundles 
(52) which are likely to be much stiffer than actin filaments. 

The case we have studied should correspond to the elastic ratchet model 
without attached filaments (13) because we have not included binding of 
filaments to the disk. However, our results appear to be at odds with the 
elastic ratchet model, which should predict a speed that depends on filament 
stiffness. One possible source of the discrepancy is that the elastic ratchet 
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model assumes that all of the force on the disk is exerted by monomers at 
the barbed end of filaments. In our simulations, roughly 40% of the force 
applied to the disk by filaments arises from monomers that are not at the 
barbed ends for the stiffer filaments we have studied. 

Why then is the speed insensitive to the bending stiffness of filaments? 
Recall that the disk repels actin, so it prefers a low concentration of actin 
near the surface. It keeps the actin concentration near its surface low by 
constantly moving forwards, away from the build-up of F-actin due to the 
action of Arp2/3. In this picture, the bending stiffness of filaments is not 
particularly important to the speed, at least at small loads. However, it is 
likely that the bending stiffness is important to other attributes of motility, 
such as the ability to withstand high loads. 

This new way of thinking about the origin of motility suggests that 
other experimental realizations of motility should be possible. Any system 
that can create a non-equilibrium, steady-state concentration profile should 
be able to develop a steady-state speed. In a real system, the mechanism is 
somewhat different because of fluid flow (53). The non-equilibrium chemical 
potential gradient resulting from the concentration gradient will lead to fluid 
flow, which will in turn push the disk. It has been understood for some 
time that a concentration gradient can lead to fluid flow which will push a 
suspended particle (53); this effect is known as "diffusiophoresis." In the case 
of actin-polymerization-driven motion of a particle such as a bacterium, bead 
or disk, the particle itself gives rise to the non-equilibrium concentration 
gradient, so the phenomenon is an example of "self-diffusiophoresis" (54). 
A recent experiment observing motility of colloids coated on one side with 
platinum that catalyzes a chemical reaction in solution is an illustration of 
a very similar phenomenon (55). 

Now that we have identified a potential mechanism, we must ask whether 
it would be significant in the real system, which we will take to be a micron- 
sized bead moving in a cell extract. The real system differs from our simula- 
tion in two very important ways. First, the parameter range is very different; 
the real system has a much lower actin concentration and depolymerization 
rate constant. Second, there is fluid flow in the real system but not in our 
simulation. To see whether the proposed mechanism is relevant to the real 
system, we estimate the speed resulting from the mechanism for realistic 
conditions within a back-of-the-envelope calculation, and compare it to the 
observed speed. If it is within an order of magnitude or so of the experimen- 
tally observed speed, our candidate is a reasonable one for the mechanism 
of motility. 

To estimate the speed, we must first estimate the concentration gradi- 
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ent of actin near the surface of the disk. Because there is a short-ranged 
repulsion between monomers and the disk, the concentration at the disk is 
approximately zero. The concentration of filaments immediately behind the 
disk, on the other hand, has been estimated by previous calculations (28, 
29) to be roughly ImM. The scale of the rise is on the order of the mesh 
size of the branched network, roughly 50 nm, so we take the concentration 
gradient to be Vc ~ lmM/50 nm. 

We must now calculate the resulting pressure gradient. A crude estimate 
is based on the ideal gas result, where Vp ~ kTVc. In the real system, this 
pressure gradient will lead to an equal and opposite pressure gradient acting 
on the fluid, which will lead to fluid flow in the actin comet tail that pushes 
the disk forwards (recall that we assume that the disk prefers to have water 
near it rather than actin). The magnitude of this flow velocity is related by 
Darcy's law to the pressure gradient of the actin via the permeability, k: 

k 

v « — Vp (2) 
V 

where i] ~ 2.4 cP is the viscosity of cell extract (56). The permeability of 
the actin comet tail can be estimated from calculations for random fiber 
networks (57) to be k ps 10~ 5 ^m 2 , but this is quite uncertain; it is only 
clear that it should be quite low. Putting this all together, we obtain a 
fluid flow speed of v ~ 1/um/s. The speed of the bead should be compa- 
rable. This speed is within an order of magnitude of the observed speed, 
which is excellent agreement despite the considerable uncertainty in the 
permeability and pressure estimates. This encouraging result suggests that 
self-diffusiophoresis is a good candidate for the origin of motility in actin- 
polymerization-driven systems. 

The proposed mechanism of motility is falsifiable by a relatively straight- 
forward experiment. According to our simulations, the key to motility lies 
in the concentration gradient of actin near the disk, which decreases as one 
approaches the disk from behind because the disk repels actin. This de- 
pends on the density of actin at the maximum, p max , which occurs roughly 
30 nm behind the surface in our simulations (see Fig. EJb)), as well as the 
density at the surface, p SU rf - I n the real system, N-WASP or Act-A at the 
surface not only activates Arp2/3 but also binds F-actin, giving rise to an 
increase in p sur f and therefore perhaps decreasing the speed. As the coverage 
of N-WASP or Act-A increases, both p max and p sur f presumably increase, 
leaving the difference relatively unaffected. This may be why the speed has 
been observed to be relatively insensitive to the coverage of Act-A (7) or 
N-WASP, at least at high coverage (34). To test our proposed mechanism, 
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we therefore propose the following experiment. Suppose one adds another 
protein to the surface, in addition to N-WASP, that binds F-actin but does 
not activate Arp2/3. By increasing the coverage of this second protein at 
fixed coverage of N-WASP, one should be able to increase p am { without af- 
fecting p m ax- This would decrease the concentration gradient, so we would 
predict that it would slow down the particle, and possibly even reverse its 
direction of motion. 

The minimal model we have presented here was designed to capture the 
most important features of the dendritic nucleation model, and we have 
tested it by reproducing results from experiments on purified proteins. One 
feature of the experimental system is missing-this is the binding of filaments 
to the protein that activates Arp2/3 complex (N-WASp, ActA, etc.). The 
next step is to incorporate specific binding of filaments to the surface. How- 
ever, we note that our success in reproducing known non-monotonic trends 
with various proteins is encouraging, and suggests that binding may not be 
essential to understanding all features of actin-based motility. 

Once we have incorporated binding, the next step will be to incorpo- 
rate bundling or crosslinking proteins and to use a curved surface. These 
extensions will allow us to study situations in which the biology has been 
perturbed, such as ActA mutants that can hop (58), bundled systems that 
still move after Arp2/3 has been removed (59), and systems that move faster 
or slower when crosslinking proteins have been added (31). 

In summary, we have conducted the first physically-consistent simula- 
tions of actin-polymerization-driven motility. These simulations are also 
the first to include semiflexible filaments and to qualitatively reproduce 
experimentally-measured, non-monotonic trends with the various proteins 
involved. Our results suggest a new picture for the mechanism of motility 
that is experimentally falsifiable. 

Methods 

Here we describe the model, which we solve numerically using Brownian 
dynamics methods. 
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Interactions 

All actin monomers, whether free or bound in filaments, are modeled as 
spheres of size a = 5nm that repel each other with a soft repulsive potential 

$R = \k Y,( R H ~ ^o) 2 , Rij <Ro = la (3) 

with K = IOO/cb^/ct 2 . Monomers within filaments interact with each other 
via a bond potential 

$s = ^^(^-^o) 2 , Rij>D = la (4) 

with K = lOOk^T / a 2 as in Eq. [3j We introduce a bending potential that 
imparts stiffness to the filament (60): 

$ B = \k b £>os(0i) - cos(0o)) 2 , (5) 

i 

We use Kb = IOO/cbT in most of our runs, but have also explored the effect 
of filament stiffness by using Kb = 01«bT and Kb = lOOOksT. In Eq. 
6 is the angle between the bond connecting monomer i — 1 to monomer i 
and the bond connecting monomer i to monomer i + 1 along a filament and 
#o = 0° (see Fig. Eta)). Note that i = 1 corresponds to the pointed end. If 
monomer i is tagged by Arp2/3 complex and is at a y-junction (37, 61, 62), 
there is also a bending potential of the same form as Eq. [5l where 9 is the 
angle between the bond connecting monomer i at the junction to monomer 
i + 1 on the branch, and the bond connecting monomer i — 1 on the parent 
filament to monomer i at the junction (Fig. IHb)). In that case, #o = 70° 
(37). 

For the moving surface, we use a flat disk of thickness a and radius 
10a (63). Monomers are repelled from the disk with a potential similar to 
Eq. [3l Note that we have not included any attractive interaction between 
filaments and the disk. As a result, the branched network is not attached 
to the disk, unlike the experimental system (29, 38-41). Model calcula- 
tions (42) suggest that the speed v (Eb) at binding energy is given by 
v(E b ) = a(E b )v(E b = 0), where a(E h ) does not depend on v(E b = 0). This 
paper focuses on the physical origin of v(E b = 0). We note that even with- 
out including binding, we are able to reproduce nontrivial, qualitative trends 
observed experimentally. Thus, it appears that binding may not be essen- 
tial to understanding all aspects of motility. Here, we have also neglected 
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crosslinking of the filaments since it is known that neither crosslinking nor 
bundling proteins are needed for motility (8). 

Branching 

In our model, the Arp2/3 complex is treated as a point particle that is gen- 
erated (activated) at the center of one side of the disk and diffuses away 
from it. We have also generated Arp2/3 at random points on one side of the 
disk and found that this makes no difference to the speed. By generating 
Arp2/3 from one side of the disk but not the other, we break symmetry. As 
a result, the branched actin network self- assembles on one side of the disk 
and drives it, on average, in a specific direction (which we define as the +z 
direction). If Arp2/3 collides with the disk it is reflected without exerting 
a force on the disk. If Arp2/3 collides with a monomer in a filament, it 
sticks to it and activates the monomer for branching. The Arp2/3 remains 
stuck to the branching monomer until the branch falls off, the branching 
monomer is depolymerized, or the Arp2/3 spontaneously dissociates. Once 
it detaches from the monomer, it is regenerated near the disk. This pro- 
cedure is designed to generate a physically reasonable Arp2/3 distribution 
near the disk surface (34, 63) without imparting forces to it as an artifact. 
We have confirmed this by running simulations with K + set to zero so that 
polymerization cannot occur. In this case, the emission of Arp2/3 from the 
disk does not lead to any motion of the disk. 

Note that we do not restrict branching to the barbed end. Because the 
number of barbed-end monomers is low compared to the total number of 
monomers in filaments, side branching (10, 37, 65) is the dominant branching 
mechanism in our model. 

Equations of motion 

All particles in our system (free monomers, monomers in filaments, Arp2/3 
and the disk) evolve according to Brownian Dynamics (Eq. [6|) (66), with 
corresponding phenomenological fricion constants C, and stochastic random 
forces F (Eq. [7|). Thus, all free monomers, filaments, Arp2/3 and the disk 
fluctuate in position due to the stochastic random forces acting on them. 
In addition, they are subjected to forces due to their interactions with each 
other: 

Ci^ = -Vi(* R + *s + *B)+-Fi (6) 
at 

(F) = 0, (Fi(t)f}(0) = 6kBT(iS(t - t')8 l3 (7) 
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The friction constant Co of an actin monomer is taken to be £o = 3irar], 
where r\ is the viscosity of the medium. The friction constant of the disk is 
taken to be £d = 2<Xo, and that for the Arp2/3 is taken to be Co- 

We convert our results to real units as follows. The unit of length in our 
model is the size of the actin monomer, a = 5nm. We take the viscosity to 
be 2.4cP, as measured experimentally for cell extracts (56). This yields a 
monomer diffusion coefficient of D = A;bT/Co = 36/xm 2 /s and a characteristic 
time unit of r = a 2 / '(2D) = 0.35/is. 

Boundary conditions 

We use periodic boundary conditions. The disk is constrained to move in 
the z-direction only. Most of our results are for a system of size 40<r x 40u x 
8O0". Unless explicitly stated otherwise, we find no discernible differences 
for system sizes 80a x 80<r x 80<r and 40<r x 40a x 160c under standard 
conditions (Table 1). 

Biochemistry 

The next step is to include the self-assembly/biochemistry of the dendritic 
nucleation model. Our algorithm for actin polymerization is similar to that 
of Gelbart et al. for nano-colloids (67). We allow polymerization only at 
the barbed end or at a branching monomer tagged by Arp2/3, and allow 
depolymerization only at the pointed end. 

Polymerization occurs when the center of a diffusing free monomer j is 
within a distance Rij of the monomer i at the growing end of a filament 
(solid rimmed circle in Fig. H{c)), such that 

(a - Sr) < Rij < a. (8) 

In addition to satisfying Eq. [8l a free monomer j must also lie within the 
angular cone 

|cos(0) -cos(0 o )| < SO (9) 

of monomer i. The angle 9 is the angle between the vector from monomer 
i — 1 preceeding monomer i on the filament to monomer i and the vector 
connecting monomer i to free monomer j (see Fig. H(c)). Here, 9q = 0°. 

Polymerization also occurs when the center of free monomer j is within 
Rij of monomer i that has been tagged by Arp2/3 as a branching monomer 
(dashed rimmed circle in Fig. Sfd)). In that case, monomer j must satisfy 
Eqs. Eland E with 9 = 70° (see Fig. 0(d)). 
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In all cases, we choose 5r and 59 in Eqs.[8]and [9] such that the potential 
energy change due to polymerization is small relative to the thermal energy 
ksT. Our choice of parameters (5r = 0.1a, 58 = 0.02) affects the effective 
polymerization rate but does not otherwise influence our results. To verify 
this, we have carried out a systematic calibration of the simulation as a 
function of polymerization rate (Appendix). 

Depolymerization is simulated using a first order rate constant, A'_. 
During each time step, At = O.OOlr, each pointed end is checked for depoly- 
merization as follows: a uniformly distributed random number between [0,1] 
is chosen and compared to the probability for dissociation during that time 
step, A_Ai. If the number is smaller than K_At the bond is broken to free 
the pointed-end monomer. Capping is treated similarly; the probability for 
a barbed end to become capped during a time step is kc+At, where kc+ 
is the pseudo-first order rate constants for capping. The probability for a 
barbed end to become uncapped is Kc-At, where Kq- is the first order 
rate for uncapping. Likewise, in each time step a branch can dissociate from 
its parent filament with probability K^At where is the debranching rate. 
Finally, we note that we do not include ADF explicitly, but instead vary the 
depolymerization rate (68). 
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Appendix: Polymerization Parameter and Calibra- 
tion 

It is well known that there is a change of free energy during polymerization; 
indeed, this is why polymerization occurs in the first place (4). The system 
gains energy by polymerizing, and gains entropy by depolymerizing. These 
free energy changes are directly related to the rate constants for polymer- 
ization and depolymerization. In the steady state system with a moving 
surface, energy is continually added to the system because the system does 
not obey detailed balance; the system is driven out of equilibrium and the 
polymerization rate is much higher relative to the depolymerization rate 
than it would be in equilibrium. In the real system, ATP hydrolysis and 
Arp2/3 activation provide this additional energy. 
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In our simulation, we do not explicitly include the free energy changes 
upon polymerization and depolymerization. Rather, we introduce rate con- 
stants that implicitly depend on those free energy changes. In our steady 
state system with a moving surface, the ratio of the polymerization rate 
to depolymerization rate has a constant value that exceeds the equilibrium 
constant, signifying that the system is out of equilibrium. 

In addition to the free energy change upon polymerization, which we take 
into account using a rate constant, we could also introduce a change in the 
mechanical energy of a filament upon polymerization. This could be done, 
for example, by storing energy in distortions of the filament. Any mechanical 
energy added through polymerization would give rise to forces that could 
lead to motility. An important question is whether polymerization in itself, 
with no mechanical energy change in the filaments due to polymerization, 
can give rise to motility. To address this question, we have designed our 
simulation model so that minimal mechanical energy is introduced into the 
filament upon polymerization. In this appendix, we will show that it is not 
necessary to introduce a mechanical energy change in the filaments in order 
to obtain motility. We note, however, that we cannot rule out the possibility 
that such a change occurs in the real system. 

It has been suggested that ATP hydrolysis may occur during Arp2/3- 
mediated polymerization at the moving surface (17, 18). In absence of direct 
experimental evidence of such a process, however, we prefer to concentrate 
on the simplest possible case, where no mechanical energy is added to the 
system even during Arp2/3-mediated polymerization, to see whether motil- 
ity and reasonable force generation can still occur. 

We cannot completely eliminate any addition of mechanical energy to 
filaments during the polymerization process. However, we can minimize it 
as follows. We have chosen the spring constant for monomer-monomer re- 
pulsion (Eq. [3j) to be the same as the spring constant holding monomers 
together in filaments (Eq. Thus, when a new bond is formed, the re- 
pulsive harmonic interaction is replaced by a full harmonic potential with 
no energy change at any value of the polymerization parameter Sr in Eq. 
However, it is impossible to avoid a mechanical energy change in the fila- 
ment due to bending of the filament (Eq. [5]) . The amount of energy change 
is determined by the parameter 59 in Eq.[9j and is nonzero as long as 59 ^ 0. 
We minimize the effect of Eq. [5] by choosing a small value for 59. To ver- 
ify that the resulting small change of bending energy does not significantly 
affect the speed, we also carry out a systematic calibration, as follows. 

The range 59 affects not only the change of bending energy stored in 
the filament due to polymerization, but also the polymerization rate itself. 
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Larger values of 59 lead to larger values of K + . Both the change in mechan- 
ical energy and the polymerization rate can, in principle, affect the speed. 
Here we check whether the dominant contribution to the change in speed 
with changing 59 arises from the change in K + , and not the change of bend- 
ing energy. To check this, we first calculate a "calibration curve" for speed 
as function of polymerization rate for 59 = 0.02. If changing 59 affects only 
the polymerization rate, then for different values of 59, corresponding to 
different polymerization rates, we should obtain speeds that lie somewhere 
on the calibration curve. 

In order to calculate the calibration curve, we must first measure the 
polymerization rate. To do this, we construct a system starting with a fixed 
small concentration of dimers, free monomers, and no disk. We turn off 
branching, capping and depolymerization and measure the rate of depletion 
of free monomers. The free monomer concentration as function of time 
is a first order decay, so we fit it to bo exp~ K+Cot where bo is the initial 
concentration of free monomers and Co is the initial concentration of dimers. 
The fitting parameter K + is the polymerization rate. The value of K + for 
our standard setup is listed in Table 1. 

The next step in calculating the calibration curve is to vary the poly- 
merization rate without changing the value of 56. This can be done without 
changing the energy of the system by introducing a probability P^ < 1 for 
capture of a monomer by the barbed end, given that the free monomer sat- 
isfies the conditions of Eqs. [8][9j In our standard runs, we use = 1, so 
we can only decrease the polymerization rate by using < 1. The result- 
ing curve for speed vs. polymerization rate is shown in Fig. [5] (rectangular 
points). 

With the calibration curve now in hand, we compute the polymeriza- 
tion rate and velocity for three different values of 59. As shown in Fig. [5j 
the speeds observed for 59 = 0.003, 0.005, and 0.015 fall on the calibration 
curve, as expected. This result demonstrates that 59 affects only the poly- 
merization rate, and that the motility is not caused by sudden changes in 
the bending energy of filaments undergoing polymerization. In other words, 
changing 59 only affects the speed through K + at small 59; there is no signif- 
icant contribution from the change of bending energy stored in the filament. 
Thus, we conclude that it is not necessary to include an explicit mechanical 
energy change upon polymerization in order to obtain motility. 

Fig. [5] shows that the speed of the disk increases linearly at low poly- 
merization rate and saturates at high polymerization rate. The satura- 
tion value of K + is related to [Arp2/3]; at low K + , the velocity is limited 
by the rate of creation of new growing ends. A straight line fit to the 
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low-polymerization-rate portion of the curve shows that a threshold poly- 
merization rate is needed to obtain a nonzero speed. This threshold rate 
yields an estimate of the critical actin concentration required for motility, 
given from Table 1. We find that the critical actin concentration is 
K^/K + ~ 2mM, consistent with what we found before in Fig. 4(d). 
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Tables 



Parameter 


zn vitro Exp. (ref) 


Simulated 


7 

bp 


5-15 u m fl3l 

. XJ i^KJ Lb 111 1 1U J 


5 //iti 


^ave 


0.1-1 urn (64, 69) 


0.1 /im 


typical bead diameter 


0.2-2 pm (56) 


0.1 /um 


viscosity (rj) 


2.4 cP (56) 


2.4 cP 


D = kBT/3-Krja 


9 — / 

36 /j,m /s 


36 /im /s 


K+ 


11.6 fzM - ^ -1 (3) 


63 pM^s' 1 


K_ 


0.3 s~l (3) 


28600 s _1 


[G-Actin] 


7 jtM (8) 


5000/iM 


if + ^G-Actin] 


270 


11 


K a 


? / uM~ i s~ i 


~K + 


K d 


0.002 s^ 1 (47) 


28600 s" 1 


[Arp2/3] 


0.1 /iM (8) 


17 


K a [Arp2/3| j. g5 -j 


N/A 


0.037 


K C + 


8 / uM~ i s~ i (65) 




K C - 


0.00042 s" 1 (65) 


s- 1 


[Cap] 


0.1 pM (8) 




/c c + = K c + [Cap] 


0.8 s - l 


14300 s" 1 



Table 1: Values of the parameters used in the simulations compared to 
those in experiments. 



27 



Figures 




Figure 1: 



Figure 1 

(a) (Color online) Disk displacement as function of time for three values 
of the bending stiffness of filaments, Kb=W00 (dotted), 100 (solid) and 
(dashed). The vertical dashed lines show the times corresponding to snap- 
shots I-III, namely 70 //s, 700 fxs and 2100 fis, respectively, after the simu- 
lation started with monomers and 5% dimers distributed randomly. Snap- 
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shots: Green spheres are monomers in filaments, red spheres are monomers 
tagged by Arp2/3 for branching, and blue spheres are the pointed-ends of 
filaments. G-actin monomers are not shown. The disk is purple. The black 
box in each frame marks the boundary of the periodic box. 

(a) 2 1 , , , 1 , , , 1 




Figure 2: 



Figure 2 

(a) Average force density on the monomers as a function of position in the 
frame of the moving disk, which is at z = 0, denoted by the vertical dashed 
line. Positive (negative) force implies monomers are being pushed to the 
right (left), (b) Total local monomer density p(z) (solid), local filament 
monomer density (dashed) and local free monomer density (dotted). 

Figure 3 

Concentration dependence of speed, (a) [Arp2/3] dependence. Solid symbols 
correspond to runs done at the standard rates shown in Table 1. Open 
symbols correspond to runs done with K_, and [G-actin] reduced by 
a factor of 10 and kc+ reduced by 5. In both cases, there is clear non- 
monotonic behavior, (b) Capping rate dependence, (c) Depolymerization 
rate dependence. Symbols correspond to the cases in which Arp2/3 protects 
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Jf_ [ms- 1 



80 - 
60 - 
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Figure 3: 



(solid) or does not protect (open) the pointed end from depolymerization. 
(d) [G-actin] dependence. 




Figure 4: 



Figure 4 

Schematic showing the definition of angles used in the bending potential 
energy in Eq. 3. (a) For a monomer i along the filament, there is a bending 
cost associated with changes of the angle 8 away from 9q = 0. (b) If monomer 
i is tagged by Arp2/3 and is at a y-j unction, there is also a bending cost 
associated with changes of the angle 6 with respect to 8q = 70°. (c) A free 
monomer j (not shown) can be added to monomer i at a growing end if its 
center is within a range Sr of separations Rij such that a— 5r < Rij < a, and 
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a range 50 of angles Oij around 6q = such that | cos 9™ — cos#o| < SO. (d) 
A free monomer j (not shown) can be added as the first monomer along a 
branch if monomer i has been tagged by Arp2/3 complex, the separation 
satisfies a — Sr < Rij < a, and the angle 6ij satisfies | cosfl^ — cos^ol < SO, 
where 9 Q = 70°. 




Figure 5: 

Figure 5 

Velocity as function of polymerization rate calibration. The data for the 
calibration points (open rectangles) are obtained with 50 = 0.02. The size 
of each rectangle corresponds to the error associated with it. 



